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We study the effects of short-time classical dynamics on 
the distribution of Coulomb blockade peak heights in a chaotic 
quantum dot. The location of one or both leads relative to 
the short unstable orbits, as well as relative to the symmetry 
lines, can have large effects on the moments and on the head 
and tail of the conductance distribution. We study these ef- 
fects analytically as a function of the stability exponent of the 
orbits involved, and also numerically using the stadium bil- 
liard as a model. The predicted behavior is robust, depending 
only on the short-time behavior of the many-body quantum 
system, and consequently insensitive to moderate-sized per- 
turbations. 



I. INTRODUCTION 

Quantum dots, semiconductor devices in which elec- 
trons are confined to live inside a two-dimensional 
mesoscopic-sized region, have generated much experi- 
mental and theoretical interest in the past decade jt). 
In the Coulomb blockade regime ||, the dot is weakly 
coupled to the outside through two narrow or tunneling 
leads, and individual resonances can be observed when 
the Fermi energy in the leads matches the energy of a 
state of N electrons in the dot. As a function of the 
Fermi energy or gate voltage one then observes a series 
of peaks in the conductance, the peak width being con- 
trolled by the temperature and the spacing between them 
by the classical charging energy required to add one more 
electron to the dot (in the experimentally typical regime 
where the level spacing is large and the intrinsic reso- 
nance width small compared with the temperature). The 
conductance peak height of the n-th resonance is then 
given by 



G n — 



where 
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and T 

am are the partial decay widths of the n-th 
resonance through each of the two leads labeled a and b. 

Each of the two partial widths is given by Fermi's 
Golden Rule as the square of a tunneling matrix ele- 
ment. This matrix element in turn is obtained (in the 



single-particle picture) by taking the overlap of the nor- 
mal derivative of the n-th dot wavefunction along the 
boundary with the electron wavefunction in the lead HQ : 



C a \(d ± V n \<t> a )\ X 



(3) 



(and similarly for r^), where 4> a is the relevant trans- 
verse wavefunction in lead a, d±^f n is the normal deriva- 
tive of the n-th dot eigenstate at the boundary, and C a is 
a constant associated with the height of the barrier (pos- 
sibly C a ^ C if the two leads are unequally coupled). 
Because tunneling will always be dominated by the lead 
mode which has the largest longitudinal energy, we may 
without loss of generality restrict ourselves to <f> a which 
is the lowest transverse energy mode of the lead [|J . For 
a smooth lead potential, this will be given by a Gaussian, 



= -(<?-9o) 2 /2o 



(4) 
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where the width a ~ Vn depends on the detailed prop- 
erties of the lead. 

Of course in reality the multi-electron state inside the 
dot is not given by a product of single-particle states, 
nor do we know the electronic Hamiltonian inside the 
dot well enough to have any realistic hope of being able 
to compute the wavefunctions ^f n . We will come back to 
these important issues in the next section. 

Many authors || have studied the behavior of the con- 
ductance peaks G n in the context of random matrix the- 
ory. There the overlaps (d± 1 $ n \4>a} are considered as ran- 
dom Gaussian variables (real or complex), and thus the 
widths T an , Tbn become y 2 random variables of one or 
two degrees of freedom in the absence or presence of a 
magnetic field, respectively. These predictions have com- 
pared favorably with the experimental data ||. In the 
present work, we extend these dynamics-free results to 
include the effects of short-time dynamics on the distri- 
bution of conductance peak heights through a ballistic 
dot. In doing so, we are following the work of Nari- 
manov, Cerruti, Baranger, and Tomsovic |0] who have 
already treated the special case of two leads placed sym- 
metrically on the horizontal-bounce orbit of a stadium 
billiard (though focusing on peak-to-peak correlations, 
rather than on the peak heights themselves, in contrast 
with the present work). Here we consider in full gener- 
ality the short-time classical effects on conductance peak 
heights, including the dependence of the peak distribu- 
tion on the stability exponents of the short orbit or orbits 
near which one or both leads may be located. We also 
disentangle the effects of symmetry lines and symmet- 
ric lead placement from the effects of short-time classical 
dynamics. 
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We note that although quantum dot experiments pro- 
vide one of the motivations for the present work, the 
effects we are considering are relevant to a wide vari- 
ety of physical situations where observable properties 
are affected by the statistical properties of wavefunc- 
tions. In particular, the conductance problem discussed 
here is formally analogous to a situation often encoun- 
tered in molecular or nuclear physics, where a reaction 
A — > B — > C takes place through an intermediate or 
transition state B living inside a metastable "well," and 
the reaction rate is then determined by the structure of 
the wavefunction B inside the well. 

The paper is organized as follows: in Section |TJ we 
briefly review some recent results concerning short-time 
dynamical effects on wavefunction intensities. We fo- 
cus in particular on the separation of scales between the 
bounce time in the dot and the time at which eigen- 
states are resolved, and on the consequent robustness of 
short-time effects on statistical wavefunction properties. 
In Section III we analyze the effects of short periodic or- 
bits (s.p.o/sj on conductance peak statistics for several 
qualitatively different situations: one lead on a s.p.o., 
two leads on the same s.p.o., and two leads on different 
s.p.o. 's for which the spectral envelopes may be in or out 
of phase with each other in the energy range of interest. 
Numerical tests of these predictions appear in Section [Tv|, 
followed by concluding remarks in Section |v|. 



II. SCARS AND WAVEFUNCTION 
INTENSITIES: BASIC RESULTS 

The scar effect is one of the most visually striking as- 
pects of quantum chaotic behavior, ft was noticed al- 
ready in the 1980's that the quantum wavefunctions of 
classically chaotic systems display anomalous enhance- 
ment and suppression of intensity in the vicinity of the 
unstable periodic orbits, contrary to the naive expecta- 
tion of wavefunction randomness and uniformity |?j . The 
early theories of this phenomenon || treated the short- 
time linearized dynamics around an unstable orbit, and 
thus made predictions about energy-smoothed spectral 
properties. More recently, theories of scars were extended 
to include long-time nonlinear recurrences, making possi- 
ble predictions about the distribution of individual eigen- 
state intensities on a given periodic orbit M. The scar 
formalism has also been adapted to study quantitatively 
wavefunction structure in systems as varied as Sinai-type 
billiards and in two-body random interaction ensembles 
in nuclear physics. 

The key result of this work that we need for the present 
analysis is that wavefunction intensities in a closed sys- 
tem are given by 



K*n|0)| 5 



r <i ! 



smooth 



(En), 



(5) 



where S^ nooth (E) is a smooth local density of states ap- 
propriate to the test state <f>, and the r n are random x 2 



variables. The smooth envelope can be determined by 
Fourier transforming the short time autocorrelation func- 
tion of 4>: 



ST ooth {E) = YT[Af mt (t)) , 



where 
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(i) = m{t))e 
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with some cutoff time To. Now as long as the cutoff Tq 
is chosen to be small compared with the Ehrenfest time, 
which scales as the inverse Lyapunov exponent of the 
system times log (kpL) (L being the size of the dot), the 
short-time dynamics of the wavepacket <f> may be deter- 
mined semiclassically, using classical motion linearized 
around the center of the wavepacket. 

In particular, let <p (which we will later identify with 
(j) a ) be a Gaussian wavepacket centered on an unstable 
periodic orbit of period P. Then we have for t = mP Jl(]] 



A 



short 



(*) = 



V cosh 0m + iQ sinh 0m 



(8) 



where > is the instability exponent for one itera- 
tion of the periodic orbit, —8 is the classical action in 
units of h (plus Maslov indices, if any), and Q is a non- 
optimality parameter. The largest recurrences (Q = 0) 
are obtained when the initial wavepacket <j> is optimally 
oriented with respect to the stable and unstable mani- 
folds of the orbit. For a lead located on a fixed periodic 
orbit, Q will in general be a function of the width a of 
the Gaussian mode in the lead (see Eq. | and Ref. @). 
More important than the analytic form of Eq. ||, how- 
ever, is the fact that for weakly unstable orbits (strictly 
<C 1, though due to numerical factors (3 « 2 is already 
in a sense "weak"), strong recurrences in A(t) persist for 
0(0 ) periods. In the energy domain, we get bumps 
in the smoothed local density of states S^ lmoth (E), of 
width 0(0) <C 1 compared to the spacing between the 
bumps, and of height 0(0^) 3> 1 compared to the 
mean. Very roughly speaking then, O(0) of all wave- 
functions are "scarred" on the periodic orbit, having in- 
tensity there O(0 _1 ) greater than the mean, while most 
of the remaining wavefunctions are "antiscarred," their 
intensity on the periodic orbit being much smaller than 
the mean. This separation of wavefunctions into scarred 
and antiscarred ones is of course only a oversimplified 
picture and one can obtain quantitatively the full dis- 
tribution of wavefunction intensities on the orbit as a 
function of the exponent 0; such a comparison between 
an analytic result for the tail of the intensity distribution 
and numerical data was in fact performed in 1 1 1 . One 
can also study the moments of the intensity distribution, 
and finds, for example, that the mean squared intensity 
approaches ir/0 times the RMT expectation for small 
(in units where the mean intensity in normalized to one). 

We emphasize the distinction between scar predictions 
and the brute-force semiclassical computation of chaotic 
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wavefunctions J12[ . In many physically interesting situa- 
tions, the Hamiltonian of the system is not known nearly 
well enough to compute individual eigenlevels and eigen- 
states either semiclassically or indeed using the full quan- 
tum machinery. What is of interest in such situations is 
not so much the detailed structure of the n-th eigenstate 
in a given sample, but rather attaining a theoretical un- 
derstanding of the statistical properties of the system. 
In a billiard (hard wall) system, changing the boundary 
by even one square wavelength far away from the pe- 
riodic orbit of interest will already destroy the detailed 
structure of individual wavefunctions on the orbit, but 
will not affect statistical properties such as the distribu- 
tion of intensities on the orbit. In the time domain we see 
this easily by recognizing that to resolve individual eigen- 
states requires following the dynamics to times of order 
of the Heisenberg time, which scales as h over the mean 
level spacing, or O(fcpL) bounce times in d = 2. Scar 
predictions, on the other hand, only require short-time 
dynamical information, on the scale of the one-bounce 
time Tg (or, to get the full effect, the classical Lyapunov 
decay time Tg / f3 if the instability exponent j3 is small) . 
Thus, deformations of the Hamiltonian will not affect 
these statistical predictions as long as the mean free path 
associated with such deformations is larger than the sam- 
ple size L (i.e. as long as we truly remain in the ballistic 
regime) . 

In the energy domain, the mean level spacing in d = 2 
goes as 1/L 2 (in units where m, % ~ 1), while the scale 
associated with scar effects (i.e. the separation between 
the bumps in the smooth local density of states enve- 
lope) is 1/Tg ~ kp/L. Thus a single bump in the LDOS 
corresponds to 0{kpL) ^> 1 level spacings, and in the 
high-energy regime (i.e. many electrons in the dot) the 
resulting wavefunction intensity statistics will be quite in- 
sensitive even to perturbations that are large compared 
with the mean level spacing. 

The preceding argument applies also to electron- 
electron scattering effects. It is known that a rather 
weak interaction beyond mean field will completely de- 
stroy level repulsion and cause the distribution of level 
spacings to approach a Gaussian form, in strong con- 
trast with the Wigner-Dyson prediction of the single- 
particle theory |13fl. The Gaussian behavior is indeed 
what is observed experimentally in such systems. This 
is not surprising because strong multi-particle effects on 
the level spacing scale require only that the interaction 
mean free path be smaller than the very large Heisen- 
berg scale ~ (kFL)TB- On the other hand, as we have 
seen, the effects on wavefunction statistics will be weak 
as long as an individual electron can freely travel across 
the device before interacting. 

It is also known (by comparing the ground state of 
a dot with TV electrons with the excited states of the 
same dot containing only N' < N electrons) that adding 
electrons to the dot changes the shape of the mean field 
potential and thus has a significant effect on the charac- 
ter of the single-particle states. Arguments very similar 



to those in the preceding paragraphs tell us that a very 
small change in the effective potential (resulting in ma- 
trix elements of the perturbation which are of order of the 
mean level spacing) is sufficient to destroy our predictive 
power for individual wavefunctions, but does not affect 
the statistical properties of these wavefunctions, which 
are associated with a much shorter time scale (the bounce 
time) and are therefore robust to any such perturbation. 
The statistical predictions would become irrelevant only 
if adding one or a few electrons to the dot changed the 
resulting potential in such a way as to completely change 
the character of the short classical trajectories. 



III. CONDUCTANCE PEAK HEIGHTS IN 
CHAOTIC SYSTEMS 

A. Generically placed leads 

Within the context of RMT, Alhas- 
sid and Lewenkopf |l4| have derived an explicit form for 
the distribution of conductance peak heights g ni allowing 
for the possibility of unequal leads (C a ^ Cb above) and 
also allowing for correlated channels. In the presence of 
time reversal symmetry,Q for two leads consisting of one 
channel each, this distribution reduces to 



P(g) 



f/r 



dr'- 



(4-7/27T 



A e -(r+r')s 



l/2». Jl/2s b y/(s a T~l/2)( Sb T' - 1/2) 




exp 



Tts*g 



(9) 



where s a and Sb are the mean partial widths through the 
two leads a and b, and in the last line we use the fact that 
the distribution depends only on the quantity s* defined 
by 



1 



1 



1 



(10) 



For equal leads, s* — s a = s&, and we also notice that in 
general the distribution of g is just the Porter-Thomas x 2 



1 For definiteness, we will consider throughout the time- 
reversal invariant situation. Of course, the calculations can 
be carried through also in the presence of a magnetic field 
(the wavefunction intensity distribution for that case has been 
studied extensively in [ fn]]ic| ]), and all predicted effects are 
qualitatively similar there. 
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distribution of one degree of freedom, with mean height 
s*/4. 

In this paper we will primarily be interested in the 
phy sical case of equal-sized leads (but see also Sec- 
tion HIE), however as we will see below the more general 



expression is a very useful starting point for studying pe- 
riodic orbit and symmetry effects. 

For two generically placed equal-sized leads, we have 
s<j = s b = s , and 



generic 



0?) 



,-2g/s 



irs g 



(11) 



The properties of this Porter-Thomas distribution are 
well-known; in particular 



so 



($)gencric — ^ cind (f/ /generic 



-ml 



(12) 



The ratio of the mean squared height to the square of 
the mean, also known as the inverse participation ratio 
(IPR) , is the simplest measure of the degree of fluctuation 
in peak intensities. 

We note also that for two leads that are symmetrically 
placed in a dot with reflection symmetry, the two partial 
widths T an and T^n are equal for each resonance n, and 
we have g n = T an /2. The peak heights are then again 
distributed according to a Porter-Thomas law, but with 
a larger mean height: 



symmetric 



(.9) 



1 



7TS0S 



-g/so 



(fj) 



•so 



symmetric 



This factor of 2 enhancement associated with perfect cor- 
relation between the two leads, and independent of any 
short-time dynamical effects, is present already in the 
case of the horizontal bounce orbit considered in Ref. || . 
We note, however, that the symmetry effect is much less 
robust than the dynamical effect, requiring the dot po- 
tential to be perfectly symmetric to better than a single 
level spacing (otherwise the even and odd states mix and 
we recover the generic result of Eq. |ll|) . 

Another independent effect, also present in the special 
example treated in [jj), is associated with the placement 
of one or both leads on the symmetry lines of the system. 
If both leads are placed on a given symmetry line, as in 
the case of the horizontal bounce orbit of the stadium, 
then only the even states produce resonances, leading to 
half as many peaks as naively expected, but with double 
the mean height of Eq. |l^: 

(ff)samc sym line = s [half expected density] . (14) 

In a system with two symmetry lines such as the Buni- 
movich stadium or the Sinai billiard, we may also con- 
sider the case of two leads on different symmetry lines. 
Then only a quarter of the eigenstates produce conduc- 
tance peaks, the two partial widths are uncorrelated but 



each is doubled with respect to the naive expectation So, 
and we obtain (cf. Eq. ^2|) 



(fj) 



diff sym lines 



so 



[quarter expected density] . (15) 



Finally, for just one lead on a symmetry line, only the 
even wavefunctions contribute with mean partial width 
2so through the symmetry line lead and sq through the 
other lead; the general expression of Eq. then leads to 



(fj) 



one on sym line 



-so 



0.343s 



(V2 + l) : 
[half expected density] 



(16) 



Thus we see that a rich diversity of conductance behavior 
may be observed simply by considering the placement of 
one or both leads with respect to the symmetry lines of 
the system. 



B. One lead on short periodic orbit 

We now put aside symmetry considerations and con- 
sider a scenario where one of the two conducting leads, 
say a, is located on a short (unstable) periodic orbit of 
the chaotic dot. Then the mean partial width through 
lead a is given by s a = s a oS snlooth (E) (see Eqs. |, |) 
in an energy range near E. Taking the two leads to be 
of equal width, s a o = S(, = sq, the effective coupling s* 
becomes energy-dependent: 



(13) S *(E) = 4s ^1 + l/y / S ,smooth (£)^ 



(17) 



Notice that the conductance is always dominated by the 
more weakly coupled lead: thus near the scarring en- 
ergies S^" 100 * 11 is strongly enhanced and the mean con- 
ductance is moderately increased (at most by a factor 
of 4) , while at the antiscarring energies S smooth is small, 
and the mean conductance may be greatly suppressed. 
Sweeping through energy, one obtains a general expres- 
sion for the distribution of conductance peak heights: 



1 



dE, 



-2g/s.(E) 



(18) 



where s*(E) is given by Eq. [j]] and S smooth (E) by the 
Fourier transform of Eq. pi Note that the energy interval 
corresponding to one oscillation of the scarring envelope 
is taken as to 27r, because we are working in units where 
h = 1 and the orbit period P is also normalized to 1. 
The distribution of Eq. |l8j may be computed numerically 
(see, for example, Fig. Hbelow) for various values of the 
stability exponent [3 (and lead non-optimality parameter 
Q, see Eq. ||). First, however, we obtain some analytic 
asymptotic results for the strongly scarred case (j3 <C 1), 
where the deviations from Porter-Thomas behavior are 
expected to be strongest. 
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The tail of the distribution P(g) will be dominated (for 
any (3) by the peak of the envelope s*(E), which coincides 
of course with the peak in S smooth (E) at E = 9 mod 2tt 
(see Eq. ||). The integral may be performed by stationary 
phase (as in Ref. Urn), to obtain 



P(g) = 



i 



-gA/2s 



2TTs g V B 



(19) 



where A = 4/s™ ax and B measures the curvature of the 
envelope at the maximum: B = 2d 2 (s~ 1 )/dE 2 . Eq. |TJ 
describes the large g behavior of the conductance peak 
distribution, where A and B are appropriate functions of 
the instability exponent (3. For small (3 one may simplify 
further and obtain A = 1 + 2^/JJC + 0{(3) and B = 
i DC -3/2p-3/2 + 0((3~ 2 ), where C w 5.24 and D « 45.1 
are numerical constants. So finally we obtain the large-g 
behavior for small (3: 



P(g) 



V2nDV 2 g ' 



-9(1+2-1 



(20) 



Notice the long tail dominated by the exponential be- 
havior exp(— g/2so) in the strongly-scarred regime, to be 
contrasted with the much shorter tail exp(— 2g/ sq) in the 
RMT ((3 -> oo) case. 

The asymptotic behavior of Eq. is plotted for in- 
stability exponent /3 = 1 as the leftmost dashed curve in 
Fig. [I], and agrees well with numerical data, represented 
by a solid curve. Thus, we see that the approximations 
used in obtaining Eq. |2C] are good already for (3 = 1, even 
though formally we have used a small-/? approximation. 
The Porter-Thomas prediction of RMT (Eq. [ll]) appears 
in Fig. [I] as a dotted line for comparison. 



one lead on periodic orbit 




two leads on different orbits 
two leads on same orbit 



FIG. 1. The tail of the conductance peak height distribu- 
tion P(g) is plotted with solid curves for (i) one lead on a 
periodic orbit with instability exponent j3 = 1, with the other 
lead generically placed; (ii) two leads on different orbits, each 
with instability exponent (3 = 1, and with in-phase spectral 
envelopes; and (iii) two leads placed on the same orbit with 
(3—1. Theoretical predictions, given by Eqs. |3f], and 
respectively, and strictly valid for g 3> ^> 1, are plotted 
as dashed curves. The mean partial width so through each 
of the two leads has been set to unity. For reference, the 
Porter-Thomas prediction of RMT (Eq. [H], to be compared 
with data sets (i) and (ii)) appears as the lower dotted curve, 
while RMT modified to include symmetry effects (Eq. [b^, to 
be compared with data set (iii)), is shown as the upper dot- 
ted curve. In the dynamics- free limit where (3 — > oo, the data 
would approach these RMT results. 

The predicted increase in the frequency of very small 
conductance peaks is even more striking. In Ref. 0] 
it was found that in between the scarring energies E = 
6 mod 2tt where the LDOS ^mooth^) ig maxnluzec i ancl 
the antiscarring energies E = 6 + n mod 2it where it is 
minimized, the LDOS at the lead follows (for sufficiently 
small instability exponent (3) the exponential law 



^smooth 



2tt 

T 



\E-e\ 



(21) 



Thus, outside of a small energy window \E— 8\ < (3\ \og(3\ 
surrounding the energy of maximal scarring, most eigen- 
states are strongly antiscarred, with mean intensity at 
the location of the lead being 5" smooth < 1. In' the 
strongly antiscarred regime where T an is small, we may 
ignore the partial width through the other lead Tb„ in 
Eq. and so the conductance g n is simply proportional 
to F Qn - S smooth {E n ). We then obtain the small-5 end 
of the conductance distribution: 



P{g) 



ir 2 g 



(22) 



which holds formally for exp(— it 2 /2f3) <j/so«/J« 1. 
The constraint exp(— n 2 /2(3) <C g/sQ is of no practical 
significance; we also note that although Eq. [2^ becomes 
exact only in the small (3 limit, quantitative agreement is 
already obtained for exponents (3 « 0.5. Eq. ^2] should be 
compared with the generic small-g behavior P ge neric(5) = 
^2/Trsog predicted by the Porter-Thomas law (Eq. [ll]) . 

The above derivation applies to an optimal lead {Q = 
in Eq. ||); for a non-optimal lead the result of Eq. 
is modified only by a Q-dependent constant, leaving the 
very distinct 1 / g scaling behavior unchanged. The small- 
g behavior of Eq. ^ is plotted (for instability exponent 
[3 = 0.2) as a dashed line in Fig. ^, and differs greatly 
from RMT expectations (the latter plotted as a dotted 
line in the same figure). 
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FIG. 2. The small-g part of the peak height distribution 
P(g) is plotted on a log-log scale for one lead on an unstable 
periodic orbit with instability exponent j3 = 0.2 (solid curve). 
The asymptotic prediction of Eq. |22| valid for g <Si (3 <Si 1 
appears as a dashed line of slope —1. For reference, the RMT 
prediction of Eq. [ll] appears as a dotted line of slope —1/2. 
Note the different power-law behavior. Again, the normaliza- 
tion is set so that the mean partial width through each lead 
is given by so = 1. The data would be almost identical for 
the cases where both leads are located on the same periodic 
orbit with = 0.2, or where they are located on two in-phase 
orbits (not shown). 

The very large fraction of small conductance peaks 
clearly must have a significant effect on the moments 
of the conductance peak distribution. From Eq. |l] we 
see that only a fraction 0(/3|log/3|) of all energies have 
^smooth -> ^ and at these energies we obtain g ~ s* ~ sq 
(see Eq. 0). At all other energies we have S smooth < 1 
and thus j ~ s, < so. All moments k > 1 of the con- 
ductance distribution then scale as 




01 log /?| 



(23) 



for sufficiently small (3. 
ments are given by 



In particular, the first two mo- 



(a) 
(g 2 ) 



0.16/3| log/3|s 
0.37/3| log P\s 2 . 



(24) 



The inverse participation ratio (IPR) is a useful dimen- 
sionless measure of the variation in heights which does 
not require one to predict theoretically the mean of the 
distribution: 



( g 2 ) 14.5 



(25) 



This result should be compared with an IPR of 3 for the 
generic Porter-Thomas distribution (see Eq. |l^). We 
see a greatly enhanced fluctuation in peak heights in the 
case where one of the leads is located in the periodic 
orbit. We also note that the numerical prefactor of « 14.5 
is valid for an optimally placed lead (Q = in Eq. 0). 



In general there is an additional prefactor which is an 
easily computable analytic function of the width of the 
Gaussian lead mode and of the monodromy matrix of the 
periodic orbit, but the important scaling behavior, i.e. 
the increase as 1//3| log /3| of the fluctuations for small (3, 
is unchanged. 




0.000001 



0.1 1 
instability exponent 



§> 0.15 




2 3 
instability exponent 

FIG. 3. The mean conductance peak height (g) is plotted 
as a function of the instability exponent /3. Three cases are 
shown: (i) one lead on an orbit of exponent f3, with the other 
lead placed generically; (ii) two leads on unrelated orbits each 
having exponent /3; and (iii) two leads on out-of phase pe- 
riodic orbits with exponent (3. For two leads on the same 
periodic orbit, the mean conductance follows the RMT pre- 
diction (with symmetry), and is independent of the exponent 
j3. (Top) Solid curves are exact; dashed curves are asymptotic 
forms valid for j3 <C 1, given by Eqs. and J37], respec- 

tively. In the dynamics-free limit f3 — > oo, all curves approach 
the RMT prediction 0.25, which is shown on the graph for 
comparison. (Bottom) Same data is shown for moderate to 
large values of the exponent (5, where the asymptotic forms 
are not applicable. 

The numerically computed mean conductance (g) for 
one lead on an orbit of instability exponent (3 appears in 
Fig. | as a function of /3 (solid curve). We observe sig- 
nificant deviations from the RMT value of 0.25 (Eq. |l2] , 
where so has been set to unity) for (3 as large as 2.0, while 
for larger values of (3 the RMT limit is approached. The 
dashed curve in Fig. ||(top) shows the asymptotic predic- 
tion of Eq. |24|, which is observed to agree well with the 
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exact results only for very small (3. 

The asymptotic behavior of Eq. ^H] appears in 
Fig. |](top) as a dashed curve, and can be compared with 
numerical data, which is plotted as a solid line. We ob- 
serve that the approximations leading to Eq. 25 do not 



lead to a quantitatively correct answer until we reach the 
very weakly unstable /3 0.1 regime. On the other hand, 
strong enhancement of the IPR compared to the RMT 
value of 3 is already clearly visible even near the mod- 
erate exponent (3 = 1, where the IPR is observed to be 
almost twice the RMT prediction. Fig. ^(bottom) shows 
the same calculation, focusing in on the IPR behavior for 
moderate values of f3, where the approximations leading 
to Eq. do not apply. 



case where the lead located on the periodic orbit also 
lies on a symmetry line of the system. As discussed in 
the previous subsection, the odd eigenstates do not then 
produce resonance peaks, while for the even states we use 
s a Q = 2so, and thus 

s,(E) = 4s ^1 + l/yj2S snw ° th (E)^j . (26) 
instead of the expression given by Eq. |l^. 



C. Leads on same periodic orbit, or on orbits related 
by symmetry 
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FIG. 4. The IPR (ratio of mean squared conductance peak 
height to the square of the mean) is plotted as a function of 
the instability exponent f) of the periodic orbit. Three cases 
are shown: (i) one lead on an orbit with exponent (3, with 
the other lead placed generically; (ii) two leads on unrelated 
periodic orbits, both with exponent f3; and (iii) two leads on 
the same periodic orbit. The RMT prediction of 3 is shown 
as a dotted line for reference; all solid curves converge to the 
RMT value in the /3 ^> 1 regime where the orbit ceases to 
be important. (Top) The solid curves are exact numerical 
results, while the dashed curves give the asymptotic /3 <C 1 
predictions of Eqs. |42| and (Bottom) Same data is 
shown for moderate values of /3, where the asymptotic forms 
are not applicable. 

The analysis in this subsection easily generalizes to the 



We now proceed to consider more generally the case 
discussed by Narimanov et al. ||, where the two leads 
are ether located on the same periodic orbit, or are lo- 
cated on orbits related by symmetry. First, we notice 
that if the two leads themselves are related by a symme- 
try, then the two wavefunction intensities and thus the 
two partial widths are identical for each resonance, just 
as in the discussion leading to Eq. [p| in the generic case. 
Less trivially, let us consider two leads located on the 
same orbit but not related by a symmetry (for example, 
one may consider the horizontal bounce in a deformed 
stadium billiard lacking a left-right symmetry, and two 
leads placed at either end of the horizontal bounce orbit). 
If the Gaussian packet corresponding to the transverse 
mode in lead b were related to the Gaussian packet asso- 
ciated with lead a by time evolution in the closed system, 
i.e. 

\<j> b ) = e-^ h \<t> a ) (27) 

for some t, then \(p a ) and \<fib) would have identical local 
densities of states, K^nl^a)! 2 = K^nl'/'o)! 2 , and thus once 
again r a „ = Tb n for each resonance n. 

More generally, even though the centers of the Gaus- 
sians \(j> a ) and |</> ) must be related by time evolution 
if they lie on the same orbit, the time-evolved version 
of \(f> a ) may have a different aspect ratio or phase space 
orientation from \4>b}- In Ref. Jill it was found, how- 
ever, that for any two optimally shaped leads, the two 
local densities of states become almost identical in the 
the limit of small instability exponent [3. Furthermore 
the minimum possible correlation between the two par- 
tial widths was shown to be 0.94 even for j3 as large as 2.3 
(corresponding to a classical stretching factor e@ = 10), 
and this correlation becomes even stronger for smaller 
j3, so that for all practical purposes in the regime where 
scarring effects are important we may take r a „ « r on for 
any two optimally oriented leads. 

In the case where either of the leads a or b are not opti- 
mally shaped so as to be aligned with the stable and un- 
stable manifolds of the periodic orbit, we have the more 
general situation 
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T an = C a (ar n + (1 - a)r an )S smooth (E n ) 
T bn = C b (ar n + (1 - a)r bn )S snwoth (E n ) , 



(28) 



where we normally consider equally-coupled leads C a — 
C b as before, and r n , r an , and r bn are independent ran- 
dom (Porter-Thomas) variables. The parameter < a < 
1 may be determined in terms of the linearized time evo- 
lution of the two wavepackets around the orbit Jl5| . In 
the weakly correlated regime a — > 0, the behavior be- 
comes identical to that of two leads on different orbits but 
having in-ph ase short-time envelopes S smooth (E). See 
Section HID below for a detailed discussion of that sit- 



uation; for now we restrict ourselves to the strongly cor- 
related case a « 1. We also note that all results valid 
for two leads on the same orbit are obviously valid also 
for leads on two orbits related by symmetry (e.g. for the 
two V-shaped orbits of the stadium). 

In the case awl, valid for two nearly optimal leads 
located on a weakly unstable orbit, the conductance peak 
height distribution reduces simply to the distribution of 
wavefunction intensities on a short periodic orbit. This 
distribution has been studied previously Jll|,0 in the 
absence of time reversal symmetry (i.e. in the presence 
of a magnetic field); here we need to develop parallel 
results for the case where a magnetic field is absent and 
the wavefunctions are therefore real. 

We begin with the moments of the distribution. At 
any given energy, the peaks heights g are distributed ac- 
cording to a Porter-Thomas distribution as in Eq. |l3| , 
but the mean height (sq/2 in the generic case) must be 
generalized to S smooth (E)s /2 (see Eq. |). The moments 
of P(g) are then given by 

{9 «) = ((S*™°^(E)) q ) ((rj-^y) , (29) 

where the average in the first factor is over energies E 
and in the second factor r n is distributed as the square 
of a Gaussian variable with variance 1. For generically lo- 
cated (though symmetrically placed) leads, S smooth (E) = 
1 independent of E and we recover the result of Eq. |l3|. 
We notice also that (S smooth (E)) = 1 by normalization, 
so the mean conductance peak height is given by 



(g) 



leads on same orbit — 



2 



(30) 



when both leads are located on the same orbit, indepen- 
dent of the stability of the orbit. This is in contrast with 
our finding in the previous subsection (Eq. ^4j) that the 
mean conductance is suppressed when only one lead is 
located on a short orbit. Furthermore, the mean squared 
intensity in this case is enhanced compared to the generic 
value: 

{(s smooth (E)) 2 ) = i^ short wi 2 

t 



cosh (3m 



7T 



(31) 



where in the first line we have used the property of 
Fourier transforms (recalling Eq. ^|), in the second line 
we have substituted from Eq. || the optimal (Q = 0) 
form of the short-time autocorrelation function, and in 
the third line we have taken the strong scarring ((3 -C 1) 
limit. Then we obtain 



(fJ 2 



3 ^> 2 



(32) 



The validity of this result is confirmed in Fig. ||(top), 
where the IPR ((g 2 ) j (g) 2 ) is plotted as a function of the 
instability exponent (3. The data for two leads on the 
same orbit appears there as the middle solid curve, which 
for (3 < 2 agrees very well with the asymptotic prediction 
of Eq. (represented by a dashed line of slope —1). The 
behavior of the same quantity for larger exponents (3 can 
be seen in Fig. ^(bottom), where we observe that already 
for (3=1 the size of peak height fluctuations is three 
times the RMT expectation, and that deviations from 
RMT are noticeable even for (3 as large as 3. 

Thus the IPR can greatly exceed the RMT value of 
3 and grows as the orbit on which the leads are placed 
becomes less unstable. This suggests that there should 
be an excess of both very large and very small peaks 
as compared with the mean (qualitatively similar to the 
result (Eq. |25| ) for only one lead on a periodic orbit). In- 
deed one can study the tail of the P(g) distribution in a 
way completely analogous to the calculation in the pre- 
vious subsection (and to the similar calculation in |ll]] 
for the non-TRS case). From Eq. 19, we obtain (us- 



From Eq. 

ing s*(E) = 2soS Smooth (E) instead of the corresponding 
form Eq. [Tt] appropriate for only one lead on a periodic 
orbit) 



P(g) = 



c 



2tt 2 D g 



P e -!3g/Cs 



(33) 



valid for small (3 and large g. As in Eq. |2(J, we see a much 
slower exponential decay in the tail than that predicted 
by RMT, but in this case the exponent is strongly (3- 
dependent and so the tail becomes ever longer as (3 — > 0. 
Eq. (for (3=1) appears in Fig. [l] as the rightmost 
dashed curve, and follows very well the exact numeri- 
cal result, which is shown as a solid curve. We see that 
the asymptotic prediction of Eq. |3^ is already very good 
even for the moderate exponent (3 = 1, while the RMT 
prediction (Eq. p"3| ) is completely wrong in the tail (com- 
pare with dotted curve), even after symmetry effects are 
included. 

This strong enhancement in the number of very large 
conductance peaks is balanced by a corresponding in- 
crease in the number of very small conductance peaks: 
the small-g part of the distribution is given again by 
Eq. ^2[ exactly as in the case of only one lead on a pe- 
riodic orbit. The small-g data for two leads on the same 
orbit is not shown here, but appears very similar to the 
case already plotted in Fig. || for only one lead on a short 
orbit. 
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The expressions obtained in this subsection need, of 
course, to be modified for situations where one of both 
leads are located on a symmetry line of the system. 
The necessary modifications are completely analogous to 
those discussed in the previous two subsections, and we 
do not go into the details here. 



D. Leads on different periodic orbits 

We now arrive at perhaps the most interesting situ- 
ation, where the two leads are located on distinct (not 
symmetry-related) periodic orbits of the classical system. 
This situation is particularly interesting because a very 
rich diversity of behavior may be obtained depending on 
the relative phase between the classical actions, in units 
of h, of the two orbits. For definiteness and for sim- 
plicity of presentation, we will focus on the case where 
two orbits have equal stability exponents (3, though the 
qualitative results obviously do not depend on this as- 
sumption. We also take the periods of the two orbits to 
be approximately equal, instead of being related by some 
other simple fraction like 1/2 or 2/3 (see the end of this 
subsection for a discussion of what happens in the case 
where the two orbit lengths are not at all simply related) . 
Then the LDOS envelopes S smooth (E) at the two leads 
are identical oscillatory functions of energy, shifted with 
respect to each other by some phase SO (see Eqs. [|, ||). 

For 59 = 0, we obtain two in-phase envelopes, and re- 
cover the a = situation of identical smooth envelopes 
with uncorrelated fluctuations under the envelopes. This 
scenario was mentioned already in the previous subsec- 
tion (Eq. ^8| ), but the discussion there was postponed 
until now. We recall that for two uncorrelated partial 
widths with the same mean, the distribution of conduc- 
tances is Porter-Thomas (Eq. [ll]), just as in the perfectly 
correlated case (Eq. |l3|), but the mean conductance in the 
uncorrelated case is only half as large as in the correlated 
case. This argument applies of course independently at 
each energy, so even after averaging over energy we obtain 
the same distribution of conductances as in the previous 
subsection, but with the mean value shifted by a factor 
of 2. To summarize, the mean is given by 



(.9) 



leads on in — phase orbits 



so 



(34) 



(cf. Eq. |3C| ), while the IPR and the small-<? part of 
the distribution (Eqs. [32], ^2|), respectively, remain un- 
changed from the previously considered case. The tail of 
the conductance peak distribution becomes 



P(g) = 



c 



2n 2 D g 



_P_ e -2t3g/Cs 



(35) 



The result of Eq. [35] is plotted for [3 — 1 as the mid- 
dle dashed curve in Fig. [|, and is observed to agree well 
with the numerical calculation, given by the correspond- 
ing solid curve. We again note that the asymptotic form 



of Eq. [35] works quantitatively even though f3 in our case 
is not small; furthermore the tail of the peak height dis- 
tribution P(g) again differs greatly from the RMT pre- 
diction (leftmost dotted curve in the same figure). 

We now consider two envelopes 5' smooth that are not in 
phase with each other. From Eq. |5| we know that each 
envelope falls off exponentially away from its peak, and 
we also know from Eq. || that when the two local den- 
sities of states are unequal, the conductance is always 
dominated by the more weakly-coupled lead (smaller 
^smooth -j rp ne max j mum conductance then occurs at en- 
ergies where S^" 100 * 11 ~ ^smooth _ rpj^g max j mum occurs 
at E = (0 a +0b)/2 mod 2ir = cf) a +5(f)/2 mod 2ir, at which 
energy S s a mooth = S s b mooth = (2?r//3) exp(-7nfy/4/?). The 
maximum possible conductance peak height is therefore 
very strongly suppressed compared with the RMT pre- 
diction: 



niax 

clout of phase orbits r "° o 



(36) 



where the out-of-phase parameter 5<j) should be taken 
between and tt. We also note that the fraction of energy 
space at which conductances this large are attained scales 
as (3 because of the exponential falloff in Eq. |l|; thus the 
mean conductance peak height scales as 



(fj) out of phase orbits ^ 



-7T50/4/3 



(37) 



suppressed by an exponential factor from the generic ex- 
pectation of Eq. |lj. The most dramatic effect arises of 
course when the two envelopes are out of phase by exactly 
tt in the energy region of interest. There we have 



fJ 



(3 



(g) 



s e 



—rr 2 /4/3 



(38) 



An exact numerical calculation for (g) is plotted as the 
lowest solid curve in Fig. |](top), and is observed to ap- 
proach the asymptotic prediction of Eq. |3^ for small /3, 
while for large (3 it of course tends to the RMT prediction 
of 0.25. Again, Fig. ||(bottom) shows in more detail the 
behavior for moderate-to-large values of (3. 

This situation described by Eq. [38] may be obtained 
in one of the following two ways. One possibility is to 
take two orbits of the same period but different Maslov 
phases, resulting in out-of phase behavior at all energies. 
Alternatively one may consider two orbits of only approx- 
imately equal period; then the two oscillating envelopes 
move in and out of phase with each other as one sweeps 
through energy by adding electrons to the dot. 

In the scenario of leads placed on two orbits with 
irrationally-related periods, we can also imagine collect- 
ing statistics over an energy window wide enough to in- 
clude both in-phase and out-of-phase behavior. We see 
in this case that the peak height will always be expo- 
nentially suppressed except at those energies where both 
LDOS envelopes S s a mooth {E) and S s b mooth (E) are near 
their respective maxima, i.e. where the wavefunction is 
scarred simultaneously at each of the two leads. Since in 
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each envelope the bumps corresponding to scarred wave- 
functions have width 0{(3) compared with the spacing 
between the bumps, and since the bumps in Sl mooth (E) 
and S b inl ° oth (E) are now assumed to be uncorrelated, we 
find that a fraction O{0 2 ) of all wavefunctions give rise 
to substantial peaks, of height O(f3~ 1 so)- Then the mean 
conductance peak height scales for small (3 as 



(<?>u 



ncorrclatcd orbits 



(3s , 



(39) 



and the higher moments of the distribution behave simi- 
larly: 



/ fe\ ffi—k k . 

\y /uncorrelated orbits H *Q ' 



' uncorrelated orbits 



/rV 



(40) 



(41) 



Again, we may measure the fluctuation in conductance 
peak heights by taking the ratio of the mean squared 
peak height to the square of the mean (IPR) . This shows 
very strong deviations from RMT expectations: 



(9) 2 



13- 



(42) 



uncorrc 



latcd orbits 



Because both leads must be scarred to obtain appreciable 
conductance, the fluctuation in conductance peak heights 
is much stronger here than in the case of only one lead 
on a periodic orbit (Eq. p5f ) or even in the case of two 
leads on the same orbit (Eq. |3^ ). 

The mean conductance peak height (g) for leads on 
two unrelated periodic orbits of instability exponent (3 
is plotted as a function of (3 as the labeled solid curve 
in Fig. Bttop) , where the (J < 1 asymptotic prediction 
of Eq. |3£j appears as a dashed line. The behavior for 
larger f3 again can be seen in Fig. ||(bottom). Similarly, 
the IPR for this case can be found plotted numerically 
by the correspondingly labeled solid curves in Fig. |](top 
and bottom); again for small j3 the result agrees with 
the asymptotic prediction of Eq. [l2] (Fig. ^(top); dashed 
line). 



E. Leads with unequal coupling 

We have been focusing throughout on the simple and 
experimentally motivated case of two equal-sized leads; 
however all of the discussion and calculations in this sec- 
tion generalize in a very straightforward way to a sce- 
nario with unequally coupled leads. Qualitatively, the 
main observation that should be added to the previous 
discussion is that the lead which is more weakly coupled 
naturally has a stronger effect on the conductance peak 
heights and their distribution, and the location of the 
more weakly coupled lead relative to classical structures 
is most important in understanding the conductance be- 
havior. 



A particularly interesting case to consider is where one 
of the leads, say 6, is coupled to the dot much more 
strongly than the other lead, so C a <C C'b in Eq. ||. Then 
Eq. |2| reduces to g n ss T an , and the conductance depends 
only on the intensit y o f the wavefunction near lead a. 
Since Sb 3> s a , Eq. [H] becomes s* = 4s a , and we ob- 
tain behavior identical to that observed previ ously for 
two symmetrically placed leads (Section [ill C ), except 
that so must be replaced everywhere by 2s a . Thus, for 
example, we obtain 



(g) = 



(43) 



(independent of the location of lead a), while the IPR 
is given by Eq. |3^ if lead a is located on an orbit of 
instability exponent (3 <C 1. The tail of the peak height 
distribution becomes 



C (3 



2tt 2 D g 



-/3g/2Cs a 



(44) 



Of course, for these results to be valid the coupling to 
lead a must be weaker than the coupling to b even near 
energies where lead a is scarred; assuming b is placed 
generically, this means C a -C /3Cb- If this condition is 
satisfied, the unequal-leads experiment may present an 
practical alternative to trying to ensure that both leads 
are centered on a periodic orbit. 



IV. NUMERICAL TESTS 

We now proceed to look at the implications of the re- 
sults of the previous section for a specific chaotic system, 
namely the Bunimovich stadium. This consists of two 
semicircular endcaps of radius 1, attached to either side 
of a rectangle of dimension 2 x 2j. The system reduces to 
a circle for 7 = 0; maximum chaos is attained at 7 = 1. 
We will focus mostly on the 7 = 1 special case, though 
scar effects are stronger at smaller values of 7 (as the 
periodic orbit instability exponents (3 become smaller). 

The (Dirichlet) wavefunctions of a billiard system such 
as the stadium can be obtained by the plane wave 
method, where at each kp one constructs the linear com- 
bination of plane waves that minimizes the integral of 
\ip 1 2 along the boundary, and then picks out those values 
of hp where this integral dips down towards zero. Of 
course, the integral never becomes strictly zero because 
of the finite number of plane waves used at any given k f ; 
however, it can easily be checked that identifying all the 
sharp minima in the boundary integral with eigenvalues 
produces the right density of states, as predicted by the 
Weyl law and higher-order corrections. For each wave- 
function obtained using this method, we can compute the 
partial widths through each of the two leads via Eq. [|, 
by simply integrating the normal derivative of the wave- 
function at the boundary, multiplied by the (Gaussian) 
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lowest transverse mode in the lead (Eq. . The conduc- 
tance peak height associated with that resonance is then 
easily obtained using Eq. |[ 

We begin by considering two leads symmetrically 
placed on the two semicircular endcaps. As we have 
seen above in Eq. [50], the mean conductance peak height 
should be independent of classical structures near the 
lead location, and this behavior is indeed observed when 
averaging over the energy range 300 < kp < 350. How- 
ever, the fluctuation in peak heights, as measured for 
example by the mean squared height, is expected to de- 
pend strongly on whether the leads are located on a short 
periodic orbit that is normally incident at the location of 
the lead (Eq. ||). 
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FIG. 5. The IPR of the conductance peak heights (mean 
squared peak height divided by the square of the mean) is 
plotted as a function of the distance of each lead from the 
center of the straight segment in a stadium billiard. The two 
leads are placed symmetrically on the two semicircular caps, 
and have width a — 1.2/vkp, which corresponds to « 3.5Af 
in the energy range considered (300 < Uf < 350). The 
amount of fluctuation is significantly enhanced when the leads 
are located on unstable periodic orbits, as compared with the 
RMT prediction of 3. The locations of the three shortest pe- 
riodic orbits having normal incidence on the semicircular end- 
caps (horizontal bounce, V-shaped orbit, and Z-shaped orbit) 
are marked on the plot. On the horizontal axis, 7 = 1.0 corre- 
sponds to the end of the straight segment, and 7 + 71-/2 « 2.57 
corresponds to the middle of the circular endcap. 

This qualitative expectation is confirmed in Fig. |^, 
where we clearly see that the mean squared peak height 
is strongly enhanced over the RMT value whenever the 
leads (of width a = 1.2/vk F « 3.4Ai?) are located on 
a short periodic orbit. Furthermore, the three shortest 
(and least unstable) orbits in the stadium having nor- 
mal incidence on the endcaps are the horizontal bounce 
(HB), the V-shaped orbit, and the Z-shaped orbit, in 
that order, and these are seen to correspond precisely 
to the three most pronounced peaks in the plot. We 
note that the periodic orbit length (as well as the insta- 
bility exponent 0) must be measured by identifying the 
four quadrants of the stadium billiard, due to symmetry: 



thus, for example, the length of the horizontal bounce 
orbit is 2 + 27 = 4. Looking at the length of the orbit 
(and instability exponent) in the full billiard, we would 
be led to underestimate the true importance of scarring 
effects. 

Notice that the leads are always located on the upper 
half of the stadium, equidistant from the center point of 
the upper straight segment. Thus, for both the horizon- 
tal bounce and the V-shaped orbit, the two leads happen 
to be on the same orbit, while in the case of the Z-shaped 
orbit, the leads would be placed on distinct orbits that 
were related by symmetry. As noted above, this differ- 
ence has no effect on the conductance behavior. 

We also see a strong increase in peak height fluctu- 
ations as the leads approach the edge of the straight 
segment. This is not surprising since the quantum be- 
havior there is strongly influenced by bouncing-ball and 
near-bouncing-ball modes |ll|. For a lead located on 
the straight segment of the boundary, the IPR will tend 
to infinity in the high-energy limit. In other words, a 
smaller and smaller fraction of all modes will give rise to 
appreciable conductance in this limit. 

To understand the conductance behavior at a quanti- 
tative level, we focus in on the most pronounced peak in 
Fig. ||, corresponding to leads located on the horizontal- 
bounce orbit, i.e. centered opposite to one another on 
the two semicircular endcaps. For 7=1, the system pre- 
viously considered in ||, the mean squared peak height 
is observed to be enhanced by a factor of 1.78 over the 
RMT prediction, as compared with an expected enhance- 
ment of 1.81 coming from Eq. |3^, where we have used 
the value f3 = 1.76 appropriate to this orbit. [We note 
that for the value of the lead width a given above, the 
lead shape is close enough to being optimal that we may 
ignore corrections associated with non-zero parameter 
Q in Eq. ||, discussed above.] Such close agreement is 
surely accidental, as one expects deviations due to the fi- 
nite number of conductance peaks being sampled, as well 
as systematic finite-energy errors associated with correc- 
tions to the semiclassical approximation and also with 
the presence of non-random bouncing ball and whisper- 
ing gallery modes. Indeed, looking at leads centered on 
the horizontal bounce orbit in the 7 = 1/2 stadium, and 
also at leads centered on the V-shaped orbit in the origi- 
nal 7 = 1 stadium, we find in each case that the observed 
mean squared peak height is about 10% lower than the 
predicted value (but still much higher than the RMT ex- 
pectation). See Table I for details. In all three cases, 
the measured value of (g 2 ) / (g) 2 is easily distinguishable 
from the IPR value of 3. 



lead locations 


(9) 


{9 2 )/{9? 




actual theory RMT 


actual 


theory 


RMT 


both on HB 


0.50 0.50 0.50 


5.33 


5.43 


3.00 


both on HB 


0.50 0.50 0.50 


6.44 


7.17 


3.00 



11 



(7 = 0.5) 






both on V 


0.50 0.50 0.50 


4.23 4.80 3.00 


HB and 
generic 


0.21 0.22 0.25 


3.70 3.73 3.00 


HB and V 


0.20 0.19 0.25 


4.11 4.50 3.00 


generic 
(symmetric) 


0.50 0.50 0.50 


3.02 3.00 3.00 



TABLE I. Mean and mean squared conductance peak 
heights are shown for several lead locations in the stadium bil- 
liard, and are compared with scar theory as well as with RMT. 
Unless stated otherwise, 7 = 1 is used for the shape of the 
billiard (see text above). We consider two cases where both 
leads are located on the same (horizontal bounce) orbit, and 
one where they are located on symmetry-related (V-shaped) 
orbits. In addition, scenarios are shown where one lead is 
located on the horizontal bounce, while the other is placed 
either on the V-shaped orbit or at a generic location (g = 2.3 
in the coordinates of Fig. ||). For reference, we also show a 
calculation where the two leads are symmetrically placed at a 
generic location; there we see that good agreement with RMT 
is obtained. 



The first two moments of P(g), as shown in Table |, 
provide useful and concise information about the distri- 
bution of conductance peak heights g, and can be easily 
used to distinguish generic lead locations from the case of 
leads located on a whost periodic orbit. More complete 
information is of course contained in the full distribution, 
as shown in Fig. ^. [The distributions there are shown as 
cumulative probabilities dg'P(g'), so as to reduce the 
effect of statistical noise.] We see, for the case of both 
leads located on the HB orbit, a long tail (Fig. ^(bottom), 
top solid curve) in good quantitative agreement with the 
prediction of scar theory (dashed curve). We also see in 
Fig. ||(top, solid curve, second from top) the increased 
number of small conductance peak heights as compared 
with RMT, and again in very good agreement with scar 
theory predictions (accompan ying d ashed curve). 

We recall also (from Section HIE) that, up to an over- 
all scaling factor of 2, exactly the same distribution would 
be obtained for only one lead located on a short orbit of 
the stadium, with the other, much more strongly coupled, 
lead located at a generic location. 
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FIG. 6. The cumulative probability J"°° dg'P(g') of having 
conductance peak height greater than g in a stadium billiard 
with two leads is plotted for several lead locations, illustrat- 
ing the various scenarios considered in the previous section: 
(i) both leads located on the horizontal bounce orbit, (ii) one 
lead on the HB orbit and the other at a generic location, (iii) 
one lead on the HB orbit and the other on the V-shaped or- 
bit (top figure only), and (iv) the two leads symmetrically 
placed at a generic location. The solid curves represent sta- 
dium data, while the dashed curves represent scar theory pre- 
dictions (Eq. [L8| etc.); for the case of generically placed leads, 
the scar theory prediction coincides with RMT (shown as a 
dotted curve). The axes are scaled so as to focus on the tail 
of the distribution in the bottom figure and the head in the 
top figure. 

We may also consider attaching two leads to the sta- 
dium so that only one o f them is located on a short orbit 
(compare Section [II B ) . We see from Table | that in this 
case the mean conductance is reduced as compared with 
RMT, while the ratio of the mean square to the square 
of the mean is enhanced, in agreement with the predic- 
tions. The full probability distribution for this case also 
appears in Fig. || (appropriately labeled solid curves, top 
and bottom), where we see a very noticeably shortened 
tail as well as an excess of small intensities, as compared 
with the previously considered case (namely, both leads 
on the HB orbit). The distribution, in the head, body, 
and tail, is in good agreement with the theoretical pre- 
diction of Eq. |l8[ 

Finally, the next-to-last line in Table | shows the con- 
ductance behavior for the case of two leads found on un- 
related short orbits, in this case the horizontal bounce 
and the V-shaped orbit. Both the suppression of the 
mean conductance and the enhancement of conductance 
fluctuations (as com pared with RMT) are observed, as 
predicted in Section [IIP . However, the size of the ef- 
fect is « 20% smaller than predicted by the theory. This 
difference may easily be ascribed to finite-energy effects 
considered above, in conjuction with statistical fluctua- 
tions. The head of the distribution is shown in Fig. ^(top, 
lowest solid curve), and follows the scar theory predic- 
tion, which again appears as a dashed curve. The tail of 
the distribution would overlap the one discussed in the 
immediately preceding paragraph (one lead on HB, with 
the other lead generic), and is not shown. 

Finally, the last line in Table |] shows the moments of 
the conductance peak height distribution for generically 
placed leads, which are seen to agree very well with RMT 
statistics. The distribution for this case is also plotted in 
Fig. ^(top and bottom), and as expected it follows the 
RMT prediction. 



V. CONCLUSION 

In this paper we have studied the statistical properties 
of transmission through a weakly open chaotic cavity. We 



have seen that RMT does a very good job in describing 
these properties when the leads are located far from clas- 
sically important structures. However, the random-wave 
hypothesis needs to be modified for a quantitatively accu- 
rate description of the transmission process when one or 
both leads are located near a classically important struc- 
ture, such as an unstable periodic orbit. The unstable 
periodic orbits have a strong effect on the distribution 
of conductance peak heights, even in the limit where the 
dwell time inside the dot becomes very large compared to 
the Lyapunov scale on which classical decay away from 
the unstable orbits takes place. This implies that short- 
time classical dynamics leaves its imprint on quantum 
behavior at arbitrarily long time scales, even though at 
these scales the classical dynamics loses all memory of its 
short time behavior. We emphasize that the effect con- 
sidered in this paper is unrelated to direct processes, and 
arises purely from quantum mechanical interference. 

We have also seen that scar-like effects on conductance 
behavior in a quantum dot are robust, and that statistical 
predictions about the distribution of peak heights can be 
made even though we do not have enough information to 
compute the individual resonances of the system, either 
quantum mechanically or semiclassically. This means in 
particular that changes in the mean field which arise as 
electrons arc added to the dot may not significantly effect 
the statistical properties of the conductance, even though 
they completely destroy all predictive power for individ- 
ual resonant wavefunctions. Similarly, electron-electron 
interactions are expected to make the single-particle pic- 
ture irrelevant for the purpose of predicting individual 
resonances or the spacing distribution, but should not af- 
fect the statistical properties of the conductance as long 
as the interaction mean free path is larger than the bal- 
listic path associated with one crossing of the dot. 

The mean conductance, higher moments, and the full 
distribution of peak heights may all be easily computed 
using scar theory ideas, for the various scenarios consid- 
ered in this paper. These include a single lead located on 
a short orbit, two leads located on distinct orbits, or two 
leads located on the same orbit. The last scenario was 
seen to be mathematically equivalent both to two leads 
located on orbits related by symmetry, and also to a situ- 
ation of unequally-coupled leads, where only the weakly 
coupled lead needs to be located on a short orbit. All 
these scenarios may be relevant to different experimental 
setups. 

Experimentally, it is of course much easier to change 
the classical dynamics inside the dot (either by adjusting 
the classical geometry or by turning on a strong mag- 
netic field) than it is to sweep through different locations 
of the lead. Thus, one may consider a situation where 
the lead locations are fixed but the classical dynamics 
is changed by adjusting either a voltage or a magnetic 
field strength. Changes in the peak height fluctuations 
are then predicted when the adjustable parameter passes 
through values which cause a periodic orbit to hit one or 
both of the leads. 
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In this paper we have not considered peak height cor- 
relations j| or the effect of finite temperature on conduc- 
tance peak heights; both of these are important effects 
which have been addressed previously by other authors. 
We have emphasized the crucial importance of including 
dynamical information for a quantitative understanding 
of conductance peak heights in a quantum dot. We also 
note that similar dynamical effects should be present in 
other physical situations where tunneling in or out of a 
chaotic well is important, and dynamical structures in- 
side the well (such as unstable periodic orbits) will influ- 
ence the distribution of coupling strengths between the 
well and the outside world. Work along these lines involv- 
ing tunneling into a smooth potential well is currently 
under development. 
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